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Abstract. Rotating shear flows, when angular momentum increases and angular 
velocity decreases as functions of radiation coordinate, are hydrodynamically stable 
under linear perturbation. The Keplerian flow is an example of such systems 
which appears in astrophysical context. Although decaying eigenmodes exhibit large 
transient energy growth of perturbation which could govern nonlinearity into the 
system, the feedback of inherent instability to generate turbulence seems questionable. 
We show that such systems exhibiting growing pseudo-eigenmodes easily reach an 
upper bound of growth rate in terms of the logarithmic norm of the involved 
nonnormal operators, thus exhibiting feedback of inherent instability. This supports 
the existence of turbulence of hydrodynamic origin in the Keplerian accretion disc 
in astrophysics. Hence, this enlightens the mismatch between the linear theory and 
experimental/observed data and helps in resolving the outstanding question of origin 
of turbulence therein. 
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1. Introduction 

Despite many efforts devoted, the origin of hydrodynamic turbulence is still poorly 
understood, because there is a significant mismatch between the predictions of linear 
theory and the experimental data (see e.g. [HE])- The mismatch between theoretical 
results and observed data is also found in astrophysical contexts, where the accretion 
flow with Keplerian angular momentum profile is a common subject. 

The hot ionized Keplerian disc flow, which is linearly stable otherwise, in presence of 
magnetic coupling is shown to exhibit magneto-rotational-instability (MRI) [3] based on 
a local analysis. Then the unstable modes were thought to lead to nonlinearity and then 
turbulence in the flow. However, it was also shown that the criterion drawn from the 
local dispersion relation is inadequate and may be misleading [H |5]. Then the question 
arises that how could a Keplerian accretion disc flow having a very low molecular 
viscosity generate turbulence, which could successively govern diffusive viscosity to 
support transfer of mass inwards and angular momentum outwards, in absence of any 
unstable mode of hydrodynamic origin? 

On the other hand, the existence of large transient growth of perturbation in plane 
Couette flow with and without the Coriolis force were demonstrated [HE]- It was 
argued that the large transient growth plausibly could lead to subcritical transition to 
turbulence in laboratory experiments and in astrophysical flows which in turn could 
lead to turbulent viscosity [6]. Note, however, that the shear flows with and without 
rotation are physically quite different. The threedimensional growth factors and then 
plausible turbulent viscosity in presence of the Coriolis force were shown to be signiflcant 
in Keplerian flow only if the elliptical vortices present in the background flow [71 [8]. 
The Coriolis effect, which absorbs the pressure fluctuation, is the main culprit to kill 
any growth of energy in the background of linear shear. Moreover, very importantly, 
the feedback of inherent instability and then sustenance of turbulence based on large 
transient growth is questionable, in lieu of linear instabilities. Therefore, the problem of 
hydrodynamic transition to turbulence in rotating shear flow, particularly when angular 
momentum increases and angular velocity decreases as functions of radial coordinate, 
e.g. in the astrophysical accretion disc, remains unsolved. While, Bech & Anderson 
[9] found turbulence persisting in numerical simulations of subcritical rotating flows for 
large enough Reynolds numbers, in general the problem of hydrodjTiamic transition to 
turbulence in rotating shear flows, particularly the astrophysical ones, has not been paid 
enough attention in the literature so far. 

Decades back Craik [10] studied the impact of rotation (including the Keplerian 
one) and hence the Coriolis effects to the flow with elliptical streamlines. Later on, 
Salhi [11] examined rotation and buoyancy effects on homogeneous shear flows and 
found similarities between rotation and stratiflcation effects. In fact, much earlier 
Bradshaw [12] discussed the formal algebraic analogy between the parameters describing 
streamline curvature and buoyancy in a turbulent flow. Very recently, the stability 
problem of unbounded rotating shear flows has been again studied in presence of uniform 
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vertical density stratification [TH] which has important astrophysical apphcations e.g. in 
geometrically thin radially stratified accretion discs. Note that indeed Johnson and 
Gammie [1^ investigated nonaxisymmetric radially stratified linearized accretion discs 
in a shearing sheet approximation neglecting vertical structure. Then they studied the 
time evolution of a plane wave perturbation co-moving with the shear fiow. The shearing 
flow with rotation and stratification has many other applications, including that in 
geophysics, as described in detail by Salhi & Cambon [13]. Interestingly, same/similar 
physics is borrowed to understand rotating shear flows in different communities with 
the use of different terminologies. While engineering community compute "rotational 
Richardson number" introduced by Bradshaw [I2] to understand linear stability of a flow 
which is proportional to the square of epicyclic frequency used by the astrophysicists 
and Rayleigh number by the fluid dynamicists in general for the same purpose. While 
positive Richardson number and Rayleigh number correspond to linear stability in 
engineering and fluid physics respectively, q < 2 (real value of epicyclic frequency) is the 
equivalent condition in astrophysics with angular velocity fl ~ r"*^ (will be discussed in 
detail in subsequent sections). 

We show that rotating flows exhibiting pseudo-eigenmodes, in absence of any 
unstable eigenmode itself, quickly reach a growth rate of the order of the logarithmic 
norm under a very small perturbation. This plausibly triggers turbulence. Note that 
the logarithmic norm, as discussed later, is an upper bound of such growth rates. 
This establishes that the intrinsic perturbation in the system is enough to sustain 
the turbulence which is a manifestation of the transient local instability. Thus, the 
turbulence of rotating shear flows and the viscosity of hydrodynamic origin in accretion 
discs is shown to be an inherent property of these systems. 

In the next section we describe hydrodynamic equations in shear flows, pseudo- 
eigenmodes, logarithmic norm and corresponding transient instability in presence of the 
Coriolis effects. Subsequently, in §3 we discuss solutions. Finally, in §4 we summarize 
with conclusions. 

2. Equations describing shear flow in presence of the Coriolis effect and the 
corresponding pseudo-eigenmode and logarithmic norm 

Let us focus on a local analysis considering a small patch of rotating shear flow whose 
unperturbed velocity proflle U corresponds to the background shear and angular velocity 
Q ^ r~^, as described earlier [2], shown schematically in Fig. [H Following earlier work 
[2\ all the variables are written in the dimensionless form. The flgure describes our choice 
of local Cartesian coordinates: x is along the radial direction, y is along the azimuthal 
or streamwise direction, and z is along the vertical direction. The unperturbed flow has 
a velocity in the y-direction and a gradient of velocity i.e. shear along the x-direction. 
The Coriolis force arised due to rotation of the patch is described by an angular velocity 
vector Q along the 2;-direction. 

For convenience, we write down the corresponding equations for evolution of 
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^^=(0,0,1) 



U=(0-x,0) 




Figure 1. Background unperturbed flow in the local comoving box. Size of arrows 
indicates the magnitude of the respective velocities. 



perturbation in terms of the set of velocity and corresponding vorticity {u, () and hence 
the Orr-Sommerfeld and Squire equations in the presence of Coriohs acceleration are 
given by 

— + f/— ^ V^n-^— + -— - — V^u (1) 

dt dy J dx"^ dy q dz Re 

— + f/—V-— V^C (2) 
dt dy J dx dz q dz Re 

where u is the x-component of velocity and ( is the corresponding component of vorticity, 
q = 3/2 for a Keplerian flow and q = 2 for a constant angular momentum flow, and 
Re is the Reynolds number. Note that in presence of linear background shear, what 
is of present interest, U = (0, —x, 0) and hence the term proportional to vanishes. 
For detailed description, see Mukhopadhyay et al. [2j. We further reduce the set of 
equations in the form of an eigenvalue equation 



where 



Q 



u 

C 

u{x, y, z, t) = u{x, t) ey.Y>[i{kyy + A;^z)], 

C{x,y,z,t) = Cix, t) exp[i{kyy + k^z)], (4) 
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and C = C{ky, kz, Re), is a differential operator. Note that because of translation 
invariance of the unperturbed flow in the y and z directions, we decompose the 
perturbations in terms of Fourier modes. 

The solution of Eqn. ([3]), describing the evolution of perturbation, can be written 
in terms of an eigenfunction expansion 

oo 

Q{x,t) = Cj exp{—iajt)Qj{x), 

*(^) - ( IS ) ■ («' 

where the complex eigenvalue of the problem cxj = a^j + icTij. Therefore, for the jth 
mode, Eqn. (|3]) reduces to 

CQj = ajQj. (6) 

To obtain the set of eigenvalues and eigenvectors, we convert the differential operator 
C into an X matrix in a finite-difference representation and then compute the 
eigenvalues and eigenvectors of the matrix. The physical parameters involved in the 
problem (mainly and the components of perturbation wave- vector ky, k^ described 
below) determine the required order N of the matrix for adequate accuracy. 

It is well known (see e.g. [2]) that all the eigenmodes are stable, i.e. ajj < for 
all j, when q < 2. However, the corresponding pseudo-eigenmodes may be unstable. 
Following previous work discussing nonrotating shear flows [Ij , let us try to understand 
it. If the linearized fluid system is driven by a signal of the form S{x,t) = e~*'^*s(x), 
then Eqn. ([3]) becomes 

^ = -tCQ + e-'^'s. (7) 

It can be verified that the response will be of the form Q{x,t) = ie~'^''^ip{x) where 
if = {al — C)~^s and / is the identity matrix. The operator (a/ — C)^-^ transforms 
the input s to the linearized flow at frequency a into the corresponding output ip. The 
maximum amplification rate that occurs in the process can be written as 

||(a/ — = maximum ^"l^^ , (8) 

where ||.{| denotes a norm of the respective quantities, in the present work 2— norm [2]. 
Now an eigenvalue of £ is a if 

Cip = criip, (9) 

where (p is the corresponding eigenvector. On the other hand, using Eqns. (|H]) and 
a can be identified as an eigenvalue of C if the perturbation with frequency a brings an 
unbounded amplification so that 

A{ky,k„R) = \\{a I - C)-^\\ = oo. (10) 
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Generalizing the above definition of an eigenvalue, we can define a as an e-pseudo- 
eigenvalue [16] of C for any e > 0, if 

A{ky,h,R) = \\iaI-C)-'\\=e-\ (11) 

when (f is the corresponding e-pseudo-eigenvector/mode such that 

||(£-a/V|| = ||6¥.||. (12) 

Naturally as e — )■ 0, and a tend to be the pure eigenmode and eigenvalue respectively. 
The set of e-pseudo-eigenvalues of an operator is called the e-pseudo-eigenspectrum. 

Therefore, we can classify the amplification in three regimes. First, if Re is equal or 
greater than a certain critical value i?ec {Re > Rec), and C has at least one eigenvalue 
in the upper half-plane of the Argand diagram, then e = and A = oo; the flow in this 
regime is linearly unstable (i.e. the system has at least one exponential growing mode). 
Second, if Re is smaller than a certain value, Reg {Re < Reg), and all the eigenvalues lie 
in the lower half-plane of the Argand diagram, then e = 1 and A = 1; this implies no 
amplification (i.e. there is no exponential growing mode). Finally, if Reg < Re < Rec 
and £ has at least one pseudo-eigenvalue in the upper half-plane (although all the 
eigenvalues still lie in the lower half-plane), then 1 > e > and 1 < A < oo; this is the 
regime of linearly stable flow with the existence of transient growth (i.e. the system has 
at least one growing pseudo-eigenmode). 

2.1. Logarithmic norm and minimum perturbation 

Now a very important related quantity is the logarithmic norm which is the upper 
bound of the possible transient growth rate. The logarithmic norm, /ip, of a square 
matrix operator A using a matrix p-norm is defined [I5] as 

r \\I + hA\\p-l 
fip := hm (13) 

for a natural number p = 1, 2, ■ ■ ■ , oo, where h is the time interval, M is the set of real 
numbers. For p = 2, it is readily shown that 

/i2(A) = A^ax (^^K^) > max {Re[X{A)]) , (14) 



2 

where A denotes an eigenvalue of the symmetric operand matrix and Re indicates real 
part. The growth of the solution Q{x, t) of Eqn. ([7]) over a non-zero time interval h can 
be upper bounded [15] as follows: 

\\Q{x,t + h)\\p<\\I + {~zC)hUQ{x,t)\\p 

+ /i|e-^"*|||s||p + C(/i2) (15) 

and the limiting maximum growth rate then can be written as 

I |Q(x, t + /.) I |p < Jim^^ (H±1z!^1^4l^ II g(^, ^ 



+ |e~'"*|||s||p + 0(/i)^ 
< fip{-tC)\\Q{x,t)\\p+ le-'^^'Wlsllp (16) 
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from which one can see that the logarithmic norm of —iC determines the upper bound of 
the transient growth of the system in Eqn. ([7]). Hence, for a suitably small perturbation 
if we obtain a pseudo-eigenvalue which is close to the logarithmic norm and is in the 
upper half plane of the Argand diagram, then the system will exhibit sustainable 
transient instability. The non-negative real parts in the pseudo-eigenspectrum of 
—iC, i.e. the values in the upper half plane of the Argand diagram for the pseudo- 
eigenspectrum of £, drive the transient growth. The pseudo-eigenvalue of —iC can be 
defined, equivalently to Eqn. ( fTTl) . as [16] 

X,{-iC) := \{-iC + E), ll^ll = e > 0,e G M. (17) 

Using the triangular inequality property of the logarithmic norm, we can write 

max [Re{X,{-iC))] < fi2{-iC + E) < fi2{-ijC) + e. (18) 

The above inequality leads to the computation of a minimum e (namely emin) for which 
the equality holds. If emin is small (compared to the matrix 2-norm of the operator, 
i.e., II — iC\\2 ) such that max [i?e(Ae(— i£))] > 0, i.e., there are pseudo-eigenvalues in 
the upper half plane of the Argand diagram, then we can conclude that the rotational 
shear flow system with the Coriolis effect exhibits transient growth sustained by small 
intrinsic perturbations. 



3. Optimal perturbations for growing pseudo-eigenmodes and logarithmic 
norms 

We solve the eigenvalue equation ([6]) with no-slip boundary conditions 
du 

M = — = C = 0, at X = ±1. (19) 
ox 

In Table 1 we enlist the optimum sets of parameters governing the perturbation that 
induce maximum growth for various values of q. We compute the maximum transient 
growth Gmax induced by the perturbation and the pseudo-eigenvalues determining the 
lower bound of the transient growth of energy: max((/m(cr)) /e)^ < Gmax, where Im 
indicates imaginary part. Interestingly, while for q = 2 the optimum perturbation giving 
rise to the maximum growth is purely vertical, for g = 1.5 it is purely two dimensional. 
For 1.5 < g < 2 it is threedimensional with ky and both nonzero. Figure [2] shows the 
optimum pseudo-eigenspectra for various values of q. Naturally for g = 2 the spectrum 
is marginally unstable and the pseudo-eigenspectra for all values of e extend to the upper 
half plane in the Argand diagram. However, as q decreases, particularly from 1.7 to 1.5, 
the pseudo-eigenvalues for particular values of e, e.g. 10~^ and 10~^'^, decrease, which 
in turn decrease the rate of maximum growth as shown in Table 1. 

As the maximum growth for g = 2 is due to the pure vertical perturbation, we 
show in Fig. [3] the pseudo-eigenspectra due to the pure vertical perturbation, giving 
rise to the maximum possible growth (maximal pure vertical perturbation), for q < 2 
as well. Corresponding parameters are given in Table 2. Interestingly, the spectra and 
thus growth for a particular q appear quite different between that for the optimum and 
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vertical perturbation (particularly for smaller values of q). This is due to the Coriolis 
effect generating epicyclic oscillations appeared in the later case, which is responsible 
for hindering growing pseudo-eigenmodes. Note that the energy growth due to pure 
vertical perturbation is bounded by the dynamics of the system and is independent of 
Re. In Fig. m we compare pseudo-eigenspectra of the flow with g = 1.5 and 1.99 for 
two different values of Re and find them very similar, reflecting their independence of 
Re under pure vertical perturbation. 

Table 3 shows the minimum perturbation (emin) for which the logarithmic norm 
rate of growth is approached. It may be noted that ^ || — '^'^Ih- Also, at g = 2 
the flow is most susceptible to transient growth since the required perturbation is the 
least. A decrease in q reduces the susceptibility to transient growth as the magnitude of 
the required perturbation increases. However, emin is still quite small even for q = 1.5, 
assuring local instability and sustained turbulence. 

4. Summary and Conclusions 

While the existence of unstable pseudo-eigenmodes favors instability and then plausible 
turbulence under finite amplitude perturbation, the question remained whether the local 
instability is sustained or not. The present analysis argues that the rotating shear flows 
with the Coriolis force not only exhibit unstable pseudo-eigenmodes under suitable 
small perturbations, but also the growth rate in the pseudo-eigenmodes approaches 
a strictly positive logarithmic norm and thus assures feedback of instability as an 
intrinsic property of the flows. Therefore, a very small amplitude of perturbation is 
able to trigger instability and then plausible turbulence in the system. In lieu of linear 
instabilities (of hydrodynamic and hydromagnetic origin) unstable pseudo-eigenmodes 
repopulating the growing disturbance assure sustained turbulence. This argues, for 
the first time, in favor of sustained hydrodynamic instability and then turbulence in 
a rotating Couette like flow, when angular momentum increases and angular velocity 
decreases with radius. This is practically the flow of many astrophysical systems, e.g. 
accretion disc around compact objects, some active galactic nuclei, star-forming systems. 
This is possible to conceive because of the non-self adjoint nature of underlying operator, 
giving rise to the non-normal eigenmodes and hence unstable pseudo-eigenspectrum 
and positive logarithmic norms. Note that molecular viscosity is too small to explain 
the observed accretion efficiencies in the systems mentioned above. Therefore, any 
viscosity is expected to arise due to turbulence in the flow. Hence, in absence of 
any linear instability in the Keplerian accretion flows, existence of small perturbations 
triggering local transient growth rates of the order of logarithmic norm, which is 
strictly positive for the associated spatial operator matrix, resolve the question of 
instability and corresponding feedback process which replenish the growth to an ultimate 
turbulence state. This is a new alternate picture of turbulence arousal in astrophysical 
flows explaining the transport of mass inward removing angular momentum outward, 
particularly when transport due to the nonlinearly developed MRI is expected to be 
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Table 1 

Parameters for optimal /maximal perturbation with Re = 2000 



Q 


ky 


kz 


Grmax 


max ([/m(cr)] /e 


1^ max [/m(cr)] 


1.5 


1.2 





21.67 


2.81 


0.165 


1.7 


1.05 


0.6 


23 


2.95 


0.159 


1.99 


0.14 


1.64 


174.43 


38.11 


4.285 X 10-2 


2 





1.66 


4661 


2199.96 


3.56 X 10-3 



Table 2 

Parameters for maximal pure vertical perturbation with Re = 2000 

q kl Gmax max {[Im{a)] /e)^ max [/m((T) 

~n IE 3^84 Tli 0.945 

1.7 2l Ol 2JL5 0.516 

1.99 L9 153.4 42.74 5.419 x lO'^ 

2 Tee 4661 2199.96 3.56 x lO'^ 

Table 3 

Minimum perturbation for logarithmic norm growth rate at Re = 2000 

q ky kz H2{.-iC,) emin II -^-^Ih 

1.5 1.2 1.1830 4.05 x 10"^ 5.1063 

~L1 1.05 0.6 1.0430 1.05 x 10'^ 5.0079 

1.99 0.14 1.64 0.3851 1.98 x lO^^ 5.0089 

2 L66 0.3632 1.86 x 10'^ 5.0087 



very small |1] in experimental setups with small magnetic Prandtl number. 
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Figure 2. Pseudo-eigenspectra of the flow under optimum perturbation when Re = 

2000 for g = 2 (uppcr-lcft), q = 1.99 (upper-right), q = 1.7 (lower-left), q = 1.5 (lower- 
right). Various curves in each figure sequentially from outside to inside respectively 
indicate results with e = 10-\ IQ-^ ■^ 10"^, IQ-^-S. 
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Figure 3. Pseudo-eigenspectra of the flow under maximal pure vertical perturbation 

when i?e = 2000 for q = 1.99 (uppcr-lcft), q = 1.7 (upper-right), q = 1.5 (lowcr-lcft). 
Various curves in each figure sequentially from outside to inside respectively indicate 
results with e = 10'^ lO'^-^, IQ-^, IQ-^-S. 
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Figure 4. Comparison of Pseudo-eigenspectra of the flows under maximal pure 

vertical perturbation between Re = 1000 (upper-left), Re = 2000 (upper-right), when 
q = 1.5; and Re = 1000 (lower-left), Re = 2000 (lower-right), when q = 1.99 Various 
curves in each flgure sequentially from outside to inside respectively indicate results 
with e = 10-\ 10-l•^ 10-2, IQ-^-S. 



